Phase diagram of the random Heisenberg antiferromagnetic spin-1 chain 
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We present a new perturbative real space renormalization group (RG) to study random quan- 
tum spin chains and other one-dimensional disordered quantum systems. The method overcomes 
problems of the original approach which fails for quantum random chains with spins larger than 
S — 1/2. Since it works even for weak disorder we are able to obtain the zero temperature phase 
diagram of the random antiferromagnetic Heisenberg spin-1 chain as a function of disorder. We find 
a random singlet phase for strong disorder and as disorder decreases, the system shows a crossover 
from a Griffiths to a disordered Haldane phase. 



The study of the effects of disorder on quantum sys- 
tems is an actual and important area of research |p]-p7|. 
Intensive work on the last decades has deepened our un- 
derstanding of the phase transitions which occur in pure 
quantum systems p8[ . Now the main effort is concen- 
trated in understanding the role of randomness in these 
transitions. This gives rise to new and interesting phe- 
nomena as the existence of Griffiths phases j|, fjj . In this 
connection random quantum spin chains have been inten- 
sive investigated. In the pure case their behavior is well 
known jljj]. Also for spin— 1/2 quantum antiferromag- 
netic chains a perturbative approach developed by Ma, 
Dasgupta and Hu (MDH) jjj and extended by Fisher 
H allows to obtain results which are essentially exact 
for this system. The picture which emerges for these 
chains is described by a random singlet phase (RSP) 
where spins are coupled in pairs over arbitrary distances. 
In the renormalization group approach this random sin- 
glet phase is governed by an infinite randomness fixed 
point ^ |7], [lTj . A straightforward extension of the MDH 
method for biquadratic spin— 1 chains has shown that in 
the Heisenberg case the perturbative RG approach may 
fail even for the case of strong disorder || . The reason is 
that in the elimination process of strong interactions, in 
which consists the MDH approach, interactions stronger 
than those eliminated are generated. This failure is bet- 
ter demonstrated when the method is extended to finite 
temperatures where it gives rise to non-physical behavior 
as negative specific heat and so on ||. Several propos- 
als have been put forward to extend the M DH method 
for quantum spin chains, with S > 1/2 jll[-|Q, without 
undisputed success. The challenge in the case of quan- 
tum integer spin chains is particularly exciting as it deals 
with the question of the fate of the Haldane phase Q 
in the presence of disorder. The existence of a gap in the 
excitation spectrum is not sufficient to guarantee the ro- 
bustness of pure chain behavior with respect to the effects 
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FIG. 1: The phase diagram of the spin— 1 random Heisenberg 
antiferromagnetic chain. Disorder is inversely related to the 
gap G of the initial rectangular distribution and Gg ~ 0.45. 



of disorder. For gapped biquadratic chains, any amount 
of disorder drives the system to a random singlet phase or 
infinite randomness fixed point . If this is not the case 
for Haldane chains there may be a unique property of in- 
teger chains which confers them a special stability with 
respect to the introduction of disorder. In fact Hida Jl(| 
using a density matrix renormalization group approach 
has not found any evidence for a RSP for spin-1 chains 
even in the presence of strong disorder, although this is 

" 13l. 



still a matter of controversy 12 



In this Letter we propose a new approach to the spin— 1 
Heisenberg antiferromagnetic chain which is an improve- 
ment of the traditional MDH perturbative renormaliza- 
tion group method. Our procedure avoids generation of 
interactions larger than those eliminated even for weak 
disorder. This provides a unique opportunity to obtain 
the phase diagram of the quantum, spin— 1 Heisenberg 
chain as a function of disorder which is shown in Fig. [l]. 
We consider here rectangular distributions of antiferro- 
magnetic interactions P(J) = (1/(1 - G))Q(1 - J)6(J- 
G). The gap G is a measure of the amount of disorder 
being inversely related to it . We find for strong disorder, 
G = 0, a random singlet phase. As G increases, i.e., dis- 
order decreases, there is a Griffiths phase characterized 
by exponents which depend on the distance to the infinite 
randomness fixed point at G = 0. This Griffiths phase 
extends up to G « 0.45 and for weak disorder there is a 
disordered Haldane phase. 



FIG. 2: The two elimination procedures as described in the 
text (Ji > J 2 ). 



The MDH method consists in finding the strongest 
interaction (Cl) between pairs of spins in the chain (see 
Fig. 2a) and treating the coupling of this pair with their 
neighbors (J\ and J2) as a perturbation. For a chain of 
spins S = 1, after elimination of the strongest coupled 
pair, the new coupling between their neighbors is given 
by, 



J' 



4 J1J2 
3 CI 



(1) 



The factor, (4/3) > 1, in this equation is the source 
of the failure of perturbation theory. Let us assume, for 
example, that the largest of the neighboring couplings 
(J\i J2) to the strongest interaction CI in the chain is J\. 
If Ji > (3/4) Cl than the new effective interaction J' is 
necessarily larger than one of those eliminated, in this 
case, than the weaker one J2. 

Our generalization of the MDH method consists in ei- 
ther of the following procedures (Fig.l). If the largest 
neighboring interaction to Cl, J\ < (3/4)f2, then we 
eliminate the strongest coupled pair obtaining an effec- 
tive interaction between the neighbors to this pair which 
is given by Eq. [l]. This effective interaction is always 
smaller than those eliminated. 

Now suppose Ji > (3/4)fi (J\ > J2). In this case, 
we consider the trio of spins S = 1 coupled by the two 
strongest interactions of the trio, J\ and Cl and solve it 
exactly (see Fig. 2b). The ground state of this trio of 
spins S = 1 is a degenerate triplet and it will be substi- 
tuted by an effective spin— 1 interacting with its neigh- 
bors through new renormalized interactions obtained by 
degenerate perturbation theory. This procedure which 
implies diagonalizing the 27X27 matrix of the trio is car- 
ried out analytically. This is important for obtaining re- 
sults on large chains and to deal with the large numbers 
of initial configurations that we use. These procedures 
guarantee that we always comply with the criterion of 
validity of perturbation theory and never, an interaction 
larger than those eliminated is generated |H| as shown 
in Fig. ||. Notice from this figure that even for the strong 
disorder case, with no gap in the original distribution of 
exchange couplings, the simple MDH procedure fails. 
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FIG. 3: Evolution of the cut-off of the random chain, as 
a function of the fraction of eliminated spins, under the two 
renormalization procedures discussed in the text. Note that 
the naive MDH process generates interaction larger than 
those eliminated, even for strong disorder (G = 0). This 
never occurs in the new procedure used here. 



We first consider the strong disorder case G — 0. This 
corresponds to the quantum critical point of the phase 
diagram where the system flows to an infinite random- 
ness fixed point. This becomes clear when we consider 
the fixed point form of the probability distribution of in- 
teractions. This is given by, 



P(J) 



a fCl 
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(2) 



The exponent a as a function of the cut-off Cl is shown 
in Fig. [|. It varies as, 



-1 
lnO' 



(3) 



This behavior characterizes the strong disorder case, G = 
0, as a random singlet phase 

Further evidence for a random singlet phase at G = 
is obtained considering the fraction of remaining active 
spins p as a function of the energy scale set by the cut-off 
Cl Q. This relation introduces a new exponent ip which 
is defined by, 



9 L |lnft|W 



(4) 



It also establishes the connection between the charac- 
teristic length L and the energy scale Cl for the case of 
logarithmic scaling. This is an extension of the usual 
definition of a dynamic exponent (f2 _1 oc r oc L z ). In 
Fig. U we show the density p = I /L as a function of the 
cut-off. From this expression we extract the exponent 
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FIG. 4: Exponent a of the fixed point, power law distribu- 
tion, Eq. |^, as a function of the cut-off Q in the low energy 
limit. 
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FIG. 5: The density of active spins as a function of the 
cut-off in the low energy limit. The inset shows the expected 
behavior for a random singlet phase with the exponent ip fn 
1/2- 



ip which takes the value ip = 1/2 characteristic of the 
random singlet phase ||. 

Finally we calculate the distribution of first gaps at 
G = 0, Hi • This is obtained starting from a given con- 
figuration of random interactions for a chain of size L and 
eliminating the spins, as described above, until a single 
pair remains. The interaction between these remaining 
spins yields the first gap A for excitation. Implement- 
ing this procedure for a large number of initial random 
configurations for chains of different sizes L yields the 
distributions P£,(logA) shown in Fig. |[ We considered 
over 10 4 initial configurations to obtain the gap distribu- 
tions. The widths of these distributions increase without 
limit as the sizes L of the chains increase, as expected for 
an infinite randomness fixed point 

According to the scaling form relating energy 
and length, Eq. ^, we expect that the distribution 
P{— log A/L^) will present a universal behavior, inde- 
pendent of the size L of the chains when plotted ver- 
sus the variable — logA/L^. This is indeed the case for 
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FIG. 6: a) Probability distribution of the first gap at the 
transition point, G = 0. The distribution become broader and 
broader with L, which signals infinite randomness behavior, 
b) Scaling plot for the gap distributions. The collapse of the 
curves is obtained for ip = 1/2 as expected for a random 
singlet phase. 



G = as shown in Fig. g for the RSP exponent ip = 1/2. 

We now decrease disorder increasing the gap in the ex- 
change distribution. In Figure |?] we show the first gap 
distributions for different degrees of disorder as charac- 
terized by the gaps G in the initial distribution of in- 
teractions. In all cases that we have investigated with 
G / 0, we find that the first gap distributions saturate 
at low energies in a form described by the expression, 
P(logA) ~ A 1 /' 2 for A — * 0. The dynamic exponent Z 
becomes independent of L for L sufficiently large. We 
have to consider large chains in order to observe this ef- 
fect. We find Zoo ~ 10.87, Z M - 1.01 and 2 m - 0.68 
for G = 0.1, G = 0.45 and G = 0.5, respectively. From 
these values of the dynamic exponent we can deduce the 
existence of a Griffiths phase extending up to Gq ~ 0.45 
where the dynamic exponent reaches the value Z = 1. 
For values of the gap G > Gq, i.e., small disorder, the 
dynamic exponent Z < 1. The distribution of first gap 
for excitations, from which low temperature thermody- 
namic properties can be deduced, implies that Z > 1 is 
required to obtain a singular behavior for these quanti- 
ties with decreasing temperature. Consequently at Go 
there is a significant change in the nature of the thermo- 
dynamic behavior of the system. The phase for G > Gg 
is a disordered Haldane phase with a pseudo-gap in the 
excitation spectrum. The thermodynamic behavior along 
the phase diagram will be explored in a future publica- 
tion ||. 

We have generalized the MDH perturbative renormal- 
ization group. This method was known to fail in the case 
of random quantum chains with spins S > 1/2 as it gen- 
erates couplings which are larger than those eliminated, 
signalling the breakdown of perturbation theory. Taking 
into account larger clusters and treating them exactly we 
were able to circumvent this problem. For the important 
case of the spin— 1 random Heisenberg antifcrromagnetic 
chain our elimination procedure gives rise to interactions 
which are always smaller than those eliminated, even for 
weak disorder. This allows us to obtain the phase dia- 
gram of this system as a function of disorder. We find 
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FIG. 7: Probability distributions of the first gap obtained 
from initial rectangular distributions of couplings with a gap 
G and different systems sizes L. For clarity not all values 
of L are shown. The solid lines represent best fits to the 
form log 10 [P(-log 10 A)] = A L - -± log 10 A. a) G = 0.1, 
Z 50 oo = 8.69, Z 6000 = 8.70, Z 7000 = 10.45, Z 8000 = 10.85 
and Z 9000 = 10.87. b) G = 0.45, Z 5000 = 0.79, Z 60 oo = 0.81, 
Z 7000 = 0.96, Z S ooo = 1.0 and Z 9000 = 1.01. c) G = 0.5, 
Zsooo = 0.55, ^6000 = 0.58, Z7000 = 0.63, Zsooo = 0.67 and 
Z9000 = 0.68. 



a random singlet phase similar to that found in spin— 1/2 
chains. For finite values of the gap G, the first gap dis- 
tribution Pi (log A) becomes, in the low energy limit, in- 
dependent of L for sufficiently large chains and is charac- 
terized by a dynamic exponent Z which depends on how 
far the system is from the infinite randomness fixed point 
at G — 0. This Griffiths phase is associated with a dy- 
namic exponent Z > 1 and is limited by the fixed point 
at G — (where Z — 00) and by Gq w 0.45, where Z 
attains the value Z = 1. For larger values of the gap, i.e., 
small disorder the system presents a disordered Haldane 
phase with a pseudogap in the spectrum for excitations. 
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for initial rectangular distributions with no gap, G = 0, 
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